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Abstract 

The doublet capture rate A1/2 of the negative muon capture in deuterium is calculated employing 
the nuclear wave functions generated from accurate nucleon-nucleon (NN) potentials constructed 
at next-to-next-to-next-to-leading order of heavy-baryon chiral perturbation theory and the weak 
meson exchange current operator derived within the same formalism. All but one of the low- 
energy constants that enter the calculation were fixed from pion-nucleon and nucleon-nucleon 
scattering data. The low-energy constant (cd), which cannot be determined from the purely 
two-nucleon data, was extracted recently from the triton /3-decay and the binding energies of the 
three- nucleon systems. The calculated values of A.1/2 show a rather large spread for the used values 
of the . Precise measurement of Ax/2 iii tbe future will not only help to constrain the value 
of d^ , but also provide a highly nontrivial test of the nuclear chiral EFT framework. Besides, 
the precise knowledge of the constant d^ will allow for consistent calculations of other two-nucleon 
weak processes, such as proton-proton fusion and solar neutrino scattering on deuterons, which are 
important for astrophysics. 

PACS numbers: 12.39.Fe; 21.45.Bc; 23.40.-s 
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I. INTRODUCTION 



The weak nuclear interaction plays crucial 
it starts the pp chain of the solar burning, 
triggered by the weak nuclear interaction [l| 

p + p — > 
p + p + e" — )■ 
p + ^He -> 
^Se + e" ^ 
«5 ^ 



le in the formation of stars in our Universe: 
In this chain, the following reactions occur, 



d + e+ + Ue, (1.1) 

d + u,, (1.2) 

'^He + e+ + z/g , (1.3) 

'Lt + V, , (1.4) 

^Be* + e+ + i/e . (1.5) 



The neutrinos produced in these reactions are messengers from the very core of the Sun, 
where the hydrogen burning occurs. Therefore, their study can provide a valuable informa- 
tion on star formation. The neutrinos, released in reaction (11. 5p have a continuous spectrum 
with the maximum energy 15 MeV and have recently been registered in the SNO detector 
[2|-01 via the reactions 

z/^ + rf z/^ + ra+p, (1.6) 
z/e + (i e~+p + p, (1.7) 

also induced by the weak nuclear interaction. As a result, the registered neutrino flux 
for the neutral current reaction fll.6p confirmed the validity of the Standard Solar Model. 
Simultaneously, the neutral current to charged current ratio [s^ established unambiguously 
the presence of an active neutrino flavor other than z/g in the observed solar neutrino flux, 
thus confirming definitely the phenomenon of the neutrino oscillations and that the neutrinos 
possess a finite mass. 

It is clear that the precise description of reactions fll.ip - fll.7l) is of fundamental value. How- 
ever, the reactions (ll.ip - fll.7p cannot be studied experimentally at present with desired 
accuracy in terrestrial conditions. In order to perceive them, one should address other weak 
processes in few-nucleon systems that are feasible in laboratories, such as 

^He + e- + Ue, (1.8) 
/i- + ^He -> + v^, (1.9) 
ix~ + d n + n+z/^. (1-10) 



Then, relying on the cosmological principle, one can apply the acquired knowledge on the 
weak nuclear interaction also to other weak processes, occurring in the extraterrestrial con- 
ditions. The reactions f ll.Sp and fll.9p have already been studied experimentally in great 
detail. The half-life of the triton is known with an accuracy ~ 0.3 %, (/Ti/2)t = (1129.6 ±3) 
s [sl, and the capture rate of muon on '^He fll.9p . Ao= 1496 ± 4 sec~^ [6, 7] is also known with 
the same accuracy. The situation with the reaction fll.lOp is less favorable so far. Indeed, 
the last measurements of the doublet capture rate provided A1/2 = 470 ± 29 sec~^ [sj and 
Ai/2 = 409 ± 40 sec~^ j^. This unfavorable situation is expected to change soon in view of 
a precision experiment planned by the MuSun Collaboration [l^ , whose goal is to measure 



2 



Ai/2 with an accuracy of ~ 1.5 %. Keeping in mind that the two-nucleon system is much 
simpler than the three- nucleon one such a precise knowledge of A1/2 will help to clarify best 
the situation in the theory of reactions triggered by weak nuclear interaction in few-nucleon 
systems. 

The main theoretical problem concerning the above quoted nuclear reactions is caused by 
the non-perturbative nature of quantum chromodynamics (QCD) at low energies. The 
way of handling this obstacle has already been outlined 50 years ago by realizing that 
the description of the electro-weak interaction with a nuclear system, containing nucleons 
and pions, should be based on the spontaneously broken global chiral symmetry SU(2)iX 



SU(2)fl, reflected in the QCD Lagrangian [ll|, |12|, see also [13|, Later on, a concept 



were 



of hidden local symmetry allowed one to extend the nuclear system to contain also heavy 
mesons The nonlinear chiral Lagrangians |17H20j served then as a starting point 

for constructing the one-boson exchange currents in the tree approximation. The reactions 
f fL6|) -f irT0D were studied with such currents in Refs. Calculations H^-il, HH, H 

performed with the nuclear wave functions derived from one-boson exchange potentials [2^ 
30|, thus providing rather consistent results. We refer to this concept as Tree Approximation 
Approach (TAA). In the more general Standard Nuclear Physics Approach (SNPA) j3l|, one 
employs realistic nucleon-nucleon (NN) potentials which usually have no relation to applied 
meson-exchange currents thus lacking the consistency. Generally, these concepts describe 
well the nuclear phenomena triggered by the electro-weak interaction up to energies ~ 1 
GeV MM. 



The crucial step allowing one to go beyond the tree approximation was made by Weinberg 
4l|. In this work, Weinberg formulated the principles of chiral perturbation theory (xPT) 
which is an effective field theory (EFT) of QCD at low energies. The effective chiral La- 
grangian is constructed by including all possible interactions between pions and nucleons 
consistent with the symmetries of QCD and, especially, the spontaneously broken approx- 



43l | allowing one to 



imate chiral symmetry. Weinberg also established counting rules j42|, 
classify the contribution of various terms of perturbative expansion in positive powers of 
g/A^, where g is a momentum or energy scale of the order of the pion mass characterizing 
a given hadronic system which is small in comparison with the chiral symmetry breaking 



scale A, 



1 GeV. 



In order to calculate reliably the capture rates and cross sections of reactions one needs to 
know accurately the current operators and the nuclear wave functions (potentials). The weak 



axial currents were studied within the EFT in Refs. |44j. |45| . At leading- (LO) and the next 



to-leading (NLO) orders, the weak nuclear current consists of the well-known single nucleon 



terms |46l . l47j . The space component of the weak axial meson-exchange current (MEG), 
which contributes to observables in all weak reactions mentioned above and is, therefore, of 
the main interest here, appears first at N^LO. It involves one unknown constant , 



d^ = di + 2d2 



1 .X 1 

- (C3 + 2C4) + - , 
6 O 



:lii) 



where the dimensionless constants q and dj are given by 



Ci = MNCi, i = l,...,A, dj 



9A 



d-n J = 1,2. 



;i.l2) 
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with Ci and dj being the low-e nerg y constants (LECs) entering the chiral Lagrangian of the 



TiNN system at the NLO |44j. |45[| . Further, M7v=0.939 GeV is the nucleon mass, gA is the 
nucleon axial vector coupling constant and denotes the pion decay constant. The linear 
combination di + 2^2 is often replaced by an effective LEG cd liF 



di + 2d2 



N 



CD 



:i.i3) 



where Av = 0.7 GeV. 



Let us note that the weak axial currents [4J, |45| were derived without referring to any 
equation describing nuclear states. The resolution of the problem of double counting in 
conjunction with the Schroedinger equation [49] provided a weak axial pion potential current, 
the presence of which ensures that the weak axial MEGs satisfy the nuclear PGAG constraint. 



q,jl^{2,q) = +([V;,p"5(l,^)] + (l ^ 2))+iUml^l{q^)M^{2,q). (1.14) 

Here V-,^ is the one-pion exchange potential, p%{l) is the one-nucleon axial charge density, 
is the pion mass, A^(g^) is the pion propagator and M"(2) is the two-nucleon pion 
absorption/production amplitude. The weak axial potential current is not a part of the 
weak axial currents of Refs. 



Eq. fll.l4p is a direct analogue of the nuclear conserved vector current constraint valid for 
the weak vector MEGs, 



g,j;(2,g-) = ([K,p"(l,g1] + (1 ^ 2)). 



;i.i5) 



The consistency of calculations requires that both potentials and currents are derived from 
the same Lagrangian. Presently, only few calculations of weak reactions performed within 
the EFT approach fulfill this requirement. Instead, one often adopts the so-called "hybrid" 
approach in which the current operator is constructed from the EFT as outlined above but 
the wave functions are generated either from the one-boson-exchange potentials of the TAA 
(30I . 50 | or from purely phenomenological potentials [5l|. Notice that all these potentials 
describe low-energy NN scattering data with x^/data ^ 1. Besides, when calculating ob- 
servables for the weak processes in the three- and more nucleon systems (see below), also 



the three-nucleon forces [52|, |53| were addressed. So far, almost all calculations aiming to 
study the weak interaction in few-nucleon systems made use of the precisely known half-life 
of the triton to extract d^ . In this way, this constant was determined in Ref. 5J] from the 
reduced Gamow- Teller (GT) matrix element for the reaction (11. 8p . and then the spectro- 
scopic factors ^pp(O) = 3.94 x (1 ± 0.004) x lO'^^ MeV b and ^hep(O) = (8.6 ± 1.3) x lO'^o 
keV b were calculated for reactions (II. ip and (II. 3p . respectively. Using the same values 
of d^ ,^the doublet capture rate A1/2 = 386 sec~^ for the reaction fll.lOp was obtained in 



56 



Ref. [55|] and the cross sections of the reactions fll.6p and (11. 7p were calculated in Ref. 
Let us note that the authors [13] employed purely phenomenological 2N AV18 51 1 and 3N 
Urbana-IX potentials. Similarly, in Ref. |57|, the capture rate Aq = 1499 ± 16 sec 



for the reaction (II. 9p was obtained based on the same potentials. In Ref. |58[, Marcucci 



and PiaruUi employed AV18 NN potential to generate the nuclear wave functions for the 
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process ^H(/i~, z/^)nn and AV18 + Urbana IX 3N potentials for the reaction ^He(/i~, t'/i)^!!. 
The weak current was taken from the xPT approach with the weak axial potential current 
of Ref. |49| added. The calculations resulted in A1/2 = 393.1 ± 0.8 s~^ and Aq = 1488 ± 9 



-1 



The problem with the above mentioned hybrid calculations is that they do not allow for a 
simultaneous determination of the LECs cd and ce which govern the short-range part of 
the 3NF at N^LO. More precisely, the constant cd enters not only the weak axial MEC, 
but also the contact and one-pion exchange part of the three-nucleon force whereas the 
constant ce controls the strength of the three-nucleon contact term jisl . 59| . These LECs 



can be determined by purely hadronic few-nucleon observables such as e.g. the triton binding 



energy and the nucleon-deuteron doublet scattering length [48| or the spectra of light nuclei 
Strong correlations between various few-nucleon observables such as e.g. the Philips 
and Tjon lines, however, often prevent from a precise determination of these LECs. On the 
other hand, a recent determination of these LECs from the triton binding energy and the 
half life |6l| seems to be more robust in this sense. 



We also mention a recent work 28|], where the values of A1/2 for the reaction (11. 101) were 
calculated within the TAA. The nuclear wave functions were derived from the Nijmegen I 
and Nijmegen 93 potentials [30| and the exchange current models were constructed within 
the approach of hidden local symmetry j20|. Taking these capture rates as input, then the 
values of were extracted in the hybrid approach, in which the nuclear wave functions 
were kept the same as in TAA and the weak axial MECs were taken from the xPT approach 
3. 



The reactions of muon capture in deuterium and in ^He have recently been studied by the 



Pisa group |62| in the SNPA- and hybrid approaches and within the xPT as well. The result- 
ing values of the capture rates are, A1/2 = (389.7 - 394.3) s~^ and Aq = (1471 - 1497) s~^. 
Subsequent calculations by this group [63] within the xPT provided A1/2 = 400±3 s^^ and 
Ao = 1494±21 s-\ We will make a more detailed comparison with these results in the next 
section. 



Recently, accurate NN potentials at N^LO by Epelbaum, Glockle and MeiBner (EGM) |6i.^ 
and Entem and Machleidt (EM) [65*1 have become available. In both cases, the parameters 
entering the calculations are standardly extracted from the fit to the NN scattering data 
and the deuteron properties, or taken from the analysis of the vrA^ scattering. 

The EGM and EM potentials differ from each other in several aspects. First, EGM adopted 
the so-called spectral function regularization [66| of the two-pion exchange contributions, 
while the analysis by EM is based upon dimensionally regularized expressions. Further 
differences can be attributed to the implementation of the momentum-space cutoff in the 
Schroedinger equation, the treatment of relativistic and isospin-breaking effects as well as 
the fitting procedure: the LECs, accompanying the contact interactions, were determined 
by EM by fitting directly to the scattering data, whereas EGM extracted them by the fit to 



the Nijmegen PWA. For more details we refer the reader to the original publications j6J, |65 



and to the review articles 1671 68 



In its turn, the simultaneous extraction of the constants cd and ce has recently been carried 



out in Ref. [6l|], where these constants have been constrained by calculations of the binding 
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energies of the three-nucleon systems and of the reduced GT matrix element of the triton 
P decay. The nuclear wave functions were generated in accurate ab initio calculations using 
both the two-nucleon N^LO EM 65| and three-nucleon N^LO force [4^, 59], whereas the 



reduced GT matrix element of the process fll.Sp has been calculated with the two-nucleon 
weak axial MECs, derived from the same xPT Lagrangian js^] as the nuclear forces. Besides, 
the calculations with two versions of the EGM potential have also been performed (without, 
however, taking into account the three-nucleon force). The resulting values of c/) cover two 
rather distinct intervals (see Table HVT) . 



In this work, we use the results of Ref. [Gjj for cd to calculate A1/2 • We apply the same 
potentials for generating the deuteron and neutron-neutron (nn) wave functions and the weak 
axial MEG operator, equivalent to that used in Ref. (6l| . thus obtaining a rather consistent 
prediction for the values of A1/2 , falling correspondingly into two rather different sets. Based 
on these results we conclude that a precise measurement of A1/2 (i) will provide a highly 
nontrivial consistency check for the chiral EFT approach and (ii) using the extracted values 
of Cd {d^ ), will allow to calculate consistently other two-nucleon weak processes, such as 
proton-proton fusion and solar neutrino scattering on deuterons, which are important for 
astrophysics. Besides it turns out [69] that enters also the capture rate for the reaction 
7r~ -|- (i — )■ 7 -|- 2 n, which is the best source of information on the nn scattering length a„„, 
see also a related work by Lensky et al. (Tof . 



In addition, we employ the weak axial pion potential current [495 ^^"^ study its influence on 
the values of A1/2 and d^ as well. As we shall see, its contribution to A1/2 is at the level 
of the weak axial exchange charge density. Therefore, it should be necessarily taken into 
account in the analysis of A1/2 • Since this current was omitted in Ref. [Glj, its influence on 
extraction of cd from the GT matrix element for the reaction (II. 8p is not clear. 

Our manuscript is organized as follows. In Section |Tll we discuss briefly the methods and 
inputs necessary for the calculations, in Section llllt we present the results and discussion 
and we conclude in Section |IVl In Appendix |Al we give a short description of our treatment 
of the nn wave functions within the K-matrix method and in Appendix |Bl the form factors, 
arising after the Fourier transformation of the weak MEGs in presence of the regulator of 
the type (12.21) . are delivered. 



II. METHODS AND INPUTS 



Here we present the necessary ingredients of the calculations of A 1/2 • They concern the 
used potentials, the method of generation of the nn wave functions, and the EFT currents. 



A. Nuclear potentials and wave functions 



In our calculations of nuclear matrix elements we start, in accordance with Ref. |6l|, by 
adopting the N^LO EM and two species of the EGM{ORC) potential with the parameters 
O = 2, R = and C = 4, 5. The meaning of these parameters is as follows: (i) O = 2 means 
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TABLE I: The scattering length a^n (in fm), and the effective range r^^ (in f^), calculated for the 
^Sq channel from the used potentials. 



^So wave 




„exp „exp 
"'nn nn 


EM 

EGM(204) 
EGM(205) 


-18.89 2.84 
-18.90 2.78 
-18.91 2.88 


-18.9 ± 0.4") 2.75 ± 0.11^) 



") Ref. [71[ ^) Ref. [7 



TABLE IL Deuteron properties derived from the EM and EGM(204) and EGM(205) chiral poten- 
tials. Here E^ is the deuteron binding energy in MeV, is the D-state probability in %, < >z 
is the deuteron magnetic moment in the nuclear Bohr magnetons, is the deuteron quadrupole 
moment in fm^, yj^r^ > is the deuteron root-mean-square radius in fm, A^ is the asymptotic 
normalization factor of the S-state in fm^^/^ and rj^ is the asymptotic D/S ratio. The magnetic 
and quadrupole moments are calculated in the impulse approximation. 

Experimental values of these quantities are: Ed=-2.224575(9)MeV [73], < /i >z=0.8574382284(94) 



7|, Q,=0.2859(3) fm^ 
%=0.0256(4) fr^. 



751], ^/<r^ >=1.9753(11) fm [76|, As=0.8846(9) fm'^s and 



Deuteron 



E, 



< f^>z Qd \/<r^ > A 



s 



Vd 



EM 

EGM(204) 
EGM(205) 



-2.224575 4.52 0.8540 0.2753 1.9750 0.8843 0.0256 
-2.218923 2.84 0.8635 0.2659 1.9856 0.8829 0.0254 
-2.223491 3.63 0.8590 0.2692 1.9754 0.8833 0.0255 



that the potential is constructed at order N^LO. (ii) The value R = provides the potential 
for the usual non-relativistic Lippmann-Schwinger equation, (iii) The value of C = 4(5) 
defines the choice of the cutoff used in the Lippmann-Schwinger equation to be 450 MeV 
(600 MeV) and of the cutoff in the spectral function representation of the two-pion potential 
to be 700 MeV (700 MeV). The nn wave functions in the continuum were generated by the 



K(R)-matrix method (50|, l78|. The basic equations are detailed in Appendix 1X1 



For a reliable evaluation of A1/2 , the potential should describe precisely the scattering length 
ann and the effective range r^n in the ^Sq channel of the nn system. In Table [U we present 
these quantities for the used potentials, together with the experimental values. We also 
present in Table HTl the properties of the used deuteron wave functions. Last but not least, 
we emphasize that in the calculations of A1/2 , the nn ^^^^Lj^ partial waves with L=0, 1,2,3 
and j/=0,l,2 are taken into account. 
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B. Weak nuclear MECs 



The weak axial MECs constructed within the HBxPT [44*] and used in our calculations 
are discussed in detail in Section 2.2 of Ref. (28|. We present in Table IIIII only the LECs 
necessary for the calculations of A1/2 . 



The LECs Cj can be obtained from vrN scattering, see Ref. (79|] and references therein. There 
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have been also attempts to determine these LECs from two-nucleon scattering data ^ ^ 

Notice that NN N^LO potentials of Refs. 64, 65[ adopt different values for the LEC C4. 
In particular, in the value of this LEC was tuned to improve the description of the 
NN data. We further emphasize that the LECs 02,3,4 receive important contributions from 
the A(1232) isobar in the process of fixing them via saturation by higher resonances, see 
Refs. HQ. 



TABLE IIL Values of the LECs a (in GeV^^) adopted in the N^LO potentials [6J,l65| 



Potential 


Cl 02 03 04 C6 


EM 
EGM 


-0.81 2.80 -3.20 5.40 3.70 
-0.81 3.28 -3.40 3.40 3.70 



In the EM model, values of C4 between 3.40 and 5.40 GeV ^ are admissible 65|. As discussed 
above, the unknown LECs di and entering the weak axial MECs, are connected with the 
LEC cd by Eq. f lLlSp . In addition, we also include the weak axial pion potential current of 
Eq. (A. 17) required by the PC AC constraint ffLTij) . 

The parameter {cd) manifests itself in the two-nucleon contact vertex of the weak axial 
current, representing a short range component of the MEC (see Fig.[l]). Besides, it is also 
present in a short-range part of the three-nucleon force. Fixing this LEC from some few- 
nucleon processes feasible in laboratory allows one to make model independent predictions 
for other weak processes (at the given order) triggered by this component of the weak axial 
current. We further emphasize that the time component of the weak axial MEC depends 
only on the known LECs C2 and C3 at the same chiral order. 

Besides these weak axial MECs, we include also the weak vector MECs of the pion range: 
the TT-pair term, the pion-in-flight term and the A excitation current of the vr range, given 
in Eqs. (A. 9), (A. 10) and (A. 11), respectively |28|. Their derivation is and properties are 
discussed in detail in Sections 2.1 and 2.2 of Ref. |28|. Let us note that the model independent 
longitudinal vr-pair and pion-in-flight terms, fixed by the low energy theorem, satisfy the 
nuclear conserved vector current constraint (11.151] , whereas the current (A. 11) is transverse 
and, therefore, model dependent. 

On the other hand, the leading order electromagnetic MECs were obtained within the 
HBxPT approach for the txN system in Ref. [83*], where they were referred to as "gener- 
alized tree graphs" . From these currents, the weak vector currents are derived by rotating 
in the isospin space. They are given in Eqs. (35) and (42) [s^. It is straightforwardly seen 
that our 7r-pair and pion-in-fiight currents coincide with the currents of Eq. (35) and the A 
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a 



b 



FIG. 1: The general structure of the two-nucleon weak axial operators; a- the long-range operator, 
b- the short-range operator. 

excitation current with the part of the current of Eq. (42) proportional to the LEG Cg, satu- 
rated by the A isobar (see Eg. (43) [83j ) *. Moreover, the same currents were obtained in the 
Appendix C of the same Ref. |83| from an effective chiral Lagrangian of the nNAptu system 
formulated in the heavy-fermion formalism. The most modern derivation of the electromag- 



netic MECs [8J, |85| , based on the chiral effective theory of the ttN system and the method 
of unitary transformation, again provides the leading order currents of the above discussed 
form (see Eq. (4.28) jsHf). These currents contain some LECs, from which so far only dig 
has been reliably fixed. The natural identification of the leading order terms, proportional 
to dis, with the standard vr-pair and pion-in-flight currents leads to 

di8 = -1/ Ami = -12.8 GeV-^. (2-1) 



On the other hand, the value of dis was fixed in Ref. [86[ by the Goldberger-Treiman dis- 
crepancy with the result (ii8=-10.14±0.45 GeV~^, which is in a good agreement with the 
value (12. ip . predicted by the low energy theorem. The other LECs, related to the transverse 
part of the discussed current are not yet well known. Natural way to extract them seems 
to be the study of pion photo- and electro-production, or capture reactions. One can also 



proceed by fixing them via saturation by the resonances, as it was done in Ref. [83 



It should be clear from our discussion that the leading order vector MECs of the pion range 
are determined by chiral invariance uniquely. Therefore, working with these MECs in this 
approximation, we can use safely our currents in the calculations done here. 



In accord with Ref. [6l|, we multiply the MECs by the regulator 

FX(g2) = e-(5/A)'" , n = 2, A = 0.5 GeV , (2.2) 



Another part of the current of Eq. (42) is proportional to the LEG Cy", which is obtained by saturation 
by the uj meson. As it is seen from Table 2 |83] , the contribution from this current to the reduced matrix 
element for the radiative neutron capture on proton is minor and we do not consider it. 
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TABLE IV: Calculated values of Aj^ , A1/2 and of the MECs effect (5MECs=Ai/2 -A1/2 (lA) (in 



s ^) for the chiral potentials EM 651] and EGM 6j]. The capture rate A1/2 differs from A1/2 by 



taking into account the weak axial pion potential current. The intervals of allowed values of cd are 
taken from Refs. 
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871 ] , the corresponding values of are derived using Eq. (jl.lip and Cj from 
Table imi 

EM"^ (EM^)- the EM potential with C4=3.4 (5.4) GeV-^ 





EM" 


EM^ 


EGM(204) 


EGM(205) 


CD 


-0.05 


0.2 


-0.3 


-0.1 


0.05 


0.33 


3.30 


4.20 




1.18 


1.44 


2.23 


2.44 


1.28 


1.57 


4.66 


5.60 


Al/2 


383.8 


387.2 


389.4 


392.4 


386.4 


389.1 


410.1 


419.1 


Al/2 


381.6 


385.1 


387.2 


390.1 


385.2 


387.9 


408.3 


417.1 


(5MECs 


4.2 


7.6 


9.8 


12.4 


2.1 


4.7 


29.9 


38.5 



keeping the currents local in configuration space. The form factors arising after the Fourier 
transform of the currents are given in Appendix [Bl 



III. RESULTS AND DISCUSSION 



The calculation of A1/2 is discussed in detail in Section 3 of Ref. [28]. Its value is obtained 
from the equation 

Al/2 = Ktat + 5A/3, (3.1) 

where h^stat and 5 A are given in Eqs. (3.2) and (3.7) [2^, respectively. The results for 
Al/2 are displayed in Table [TVl As it is seen from Table HVt the values of A1/2 are located 
in two rather distinct intervals. For the NN potentials EM and EGM(204), they are in the 
interval 383.8 s~^ < A1/2 < 392.4 s~^, thus featuring a spread of about 2%. The spread of 
the values of A1/2 obtained in the EGM(205) model is about the same, but these values are 
larger by ^ 6 % and the MECs effect is also much larger. If one would like to get within 
the model EGM(205) the value A1/2 =389.1 s~^, one should take =2.35 corresponding to 
cd=1.08 which is far out of the interval 3.30 < cd < 4.20 allowed by the reduced GT matrix 
element of the triton /3 decay. 

One can see from Table IIVI that the presence of the weak axial pion potential current 
suppresses A1/2 up to 0.5 %. E.g., in the case of the model EM'' and for the value {cd) 
=2.23 (-0.3), one gets A1/2 =389.4 s~^, whereas A1/2 =387.2 s"^. This difference in the 
capture rate can be reflected in the change of d^ [cd) = 2.07 (-0.45), which means the 
change of ~ 7.5 % (44 %) in d^ (c^). 

Let us now discuss for the same model EM'' the variation of extracted d^ (cd), if the sup- 
posed experimental value of A1/2 =390 ±6 s^^ is measured with an accuracy of 1.5 %. The 
extracted value of d^ =2.27±0.44, to which corresponds Ci:)=-0.26±0.42. It is seen that the 
value of d^ (co) would be known with an accuracy of 20 % (140 %). With the weak axial 
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pion potential current included, the interval of (co) would be shifted to the negative 
values by ~ 7.5 % (44 %) as a whole. 

In Table |Vl we give individual contributions of various components of the currents, the 
multipoles of which can be found in Appendix C of Ref. [28[. For the EGM version, the 
contribution of the weak axial pion potential current X4(p.p.c.) is about the same as that 
of the time component of the weak axial MECs, and for the EM potentials it is twice as 
large. Notice that in the calculations of [28| based on the Nijmegen I potential js^l, the 
suppression of A1/2 by ~ 1 % was obtained for this current, whereas in Ref. [ssf no effect has 
been found based on the nuclear wave functions generated by the Argonne fig potential (sij . 
For the EGM (205) version, the contribution of the short-range weak axial MECs, Ja(EFT), 
is strongly enhanced as a consequence of the large value of . 

In Table IVIt we compare the results of our lA calculations for the model EM^ with the 



second row of Table II of Ref. [63[. As it is seen, the agreement between the calculations is 
very good. On the other hand, the comparison of the final result A1/2 =392.0±2.3 s~^ [62| 
with our Table HVl shows reasonable agreement in calculations for the potential model EM^, 
whereas the result A1/2 =399±3 s~^ 63|| indicates the difference of ~ 2 - 3 %. In general, 



our results exhibit much more model dependence of A1/2 , than admitted in Refs. |62|, |63 



Let us note that our results for A1/2 contain the factor 1.028(4) taking into account the 
radiative corrections. These corrections have not yet been calculated within the EFT. The 
result of "classical" calculations j47|] for the inner correction is 0.024(4), valid for the muon 
capture in any nucleus. The vacuum polarization correction to the muon bound state wave 
function must be evaluated for every nucleus separately. For the muon capture in deuteron, 
this correction is 1.80 a;/7r=0. 0042 [88]. Then the overall radiative corrections for the muon 
capture in deuteron turn out to be the same as for the muon capture on proton. 



TABLE V: Contributions of various components of the currents to A1/2 (in s ^) obtained using 
chiral potentials EM 65| and EGM 64| for certain values of . lA(^S'o)- the contribution to 



Ai/2 from the nn final state and the one-body currents; -|- j^- with the time component of the 



weak axial MECs added; -|-X4(EFT)- with the space component of the EFT weak axial MECs [44 1 
added; +jv- with the space component of the weak vector MECs added; lA(full)- the contribution 
to Ai/2 from all the nn partial waves taken into account and the one-body currents; -|- ja- with 
the weak axial MECs 4J] added; +jv- with the space component of the weak vector MECs added; 
+Jvi(p-P-c.) - with the weak axial pion potential current added. The values of the doublet capture 
rate in the last column present in fact A1/2 ■ 



Potential 




IA(i5o) +fA +X4(EFT) +JV 


lA(full) +JA +JV 


+Ja(p-p-c.) 


EGM(204) 
EGM(205) 
EM"^ 
EM'' 


1.57 
5.13 
1.18 
2.23 


245.2 244.1 248.3 249.8 

241.8 240.7 274.7 276.6 

240.9 240.0 243.0 244.9 
240.9 240.0 248.5 250.5 


384.3 387.3 389.1 

379.5 412.4 414.6 

379.6 381.6 383.8 
379.6 387.2 389.4 


387.9 
412.7 
381.6 
387.7 
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TABLE VI: Contributions of the nn partial waves to (in s calculated for the model EM'' 
and with the lA currents. The third row contains the second row of Table II [63]. 



This work 
Ref. [631 



5K^/2{^So) Mi/2('Po) Ml/2('^l) 5^1/2{^P2) (^Al/2(^^2) 5A^/2{^F2) Ay2 



240.9 22.0 40.9 69.4 6.0 0.4 379.6 

238.8 21.1 44.0 72.4 4.5 0.9 381.7 



IV. CONCLUSIONS 



In this work, we studied the doublet capture rate, A1/2 , for the reaction fi~{d, 2n)z/^ f ll.lOp 
within the framework of chiral EFT. The nuclear wave functions of the initial and final 
states were generated from the N^LO chiral potentials of Entem and Machleidt [65[ and of 



Epelbaum, Glockle and MeiBner 6J|. The properties of the ^Sq nn state and of the deuteron 
are presented in Table |T] and Table |IIl respectively. The employed weak axial MECs were 
also constructed within the same EFT scheme The corresponding LECs, taken 



consistently with the potentials, are listed in Table IIIII The short-range part of the weak 
axial MECs contains one unknown LEG, [cd), which has been constrained earlier j6l| by 
Gazit, Quaglioni and Navratil in the study of the GT matrix element of the triton P decay 
and of the static properties of the 3N system, using the same NN potentials and weak axial 
MEGs. With the constant , constrained in this way and with the same regulator (12. 2p 
entering the weak MEGs, we performed rather consistent calculations of A1/2 • As it is seen 
from Table IIVI the results are located in two rather distinct intervals, in which the spread 
of Ai/2 is about 2 %. It follows from this result that (i) studying only the triton /3 decay or 
the muon capture in deuterium might be insufficient for a precise determination of d^ , (ii) 
combined investigation of both reactions can not only provide precise value of d^ , but will 
allow for a highly nontrivial test of the chiral EFT framework. 

In addition, we studied the influence of the weak axial pion potential current on the A1/2 and 
d^ . This current is not a part of the weak axial MEGs derived within the HBxPT in 



Refs. |44j . |45|, and was not taken into account in the calculations |61| either. As it is seen 
from our Table IIVI and Table IVj its contribution to A1/2 is non-negligible and can change 
the value of d^ (cd) up to 7.5 % (44 %). In our opinion, the full consistency of calculations 
cannot be achieved without including this current also in the study of the reaction (11. 8p . 

What are the possible origins of the observed fairly large spread in the calculated values of 
d^ when using the EGM NN potential of Ref. [6j]? First of all, the value of the LEG cd was 
obtained in Ref. {6l| without taking into account the 3NF in the case of the EGM potentials. 
Given the significant underbinding of the triton for N^LO NN potentials of Ref. [6^ without 
the inclusion of the 3NF ^, the obtained values for cd in the case of the EGM NN potentials 
should be taken with care. Furthermore, while the applied regulator keeps the calculations 
in the configuration space local, the regulator adopted in the N'^LO NN potentials [g^, 65| 



^ The underbinding is of the order of 1 ... 1.5 MeV which is larger than in the case of the N^LO potentials. 
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Fa(p',p) = e-''"/^'-^'/^', A = 0.5GeV, (4.1) 

where p [p') is the initial (final) nucleon momentum in the center-of-mass system. Clearly, 
this regulator leads to strong non-localities in configuration space. Strictly speaking, only 
calculations with this type of the regulator can be considered as fully consistent. It is 
desirable to carry out such calculations in the future in order to see how the form of the 
regulator affects the results for A1/2 and the values of extracted constant (cd)- Another 
possible source of inconsistency emerges from using two- and three-nucleon forces and the 
exchange currents at different orders in the chiral expansion. The N^LO NN potentials 
should be accompanied with the 3NFs calculated at the same order in the chiral expansion. 
Only recently, the N^LO corrections to the 3NF have become available 89|- Moreover, 



the results of the pioneering work [92| indicate that the inclusion of the long-range N^LO 
corrections in the 3NF might have strong impact on the values of the LECs cd and ce- This 
issue needs to be investigated in the future. Finally, we applied the weak vector MECs in the 



lea ding order only. Evidently, the estimation of higher order terms, derived in Refs. [83 



85l . l93l |94| |. is mandatory in future calculations. Last but not least, it would be interesting 
to carry out the calculations using the N^LO NN potentials and the corresponding 3NF. 
This would help to obtain a reliable estimation of the theoretical uncertainty for A1/2 ■ 
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Appendix A: The nn wave functions 



Here we discuss in more detail the derivation of the nn wave functions in the formahsm 
of the K-matrix, referring essentially to Ch. 6 and Ch. 7 of Ref. 78|] and to Appendix A of 
Ref. I5C 



The basic equation for the nn scattering wave function in terms of the T-matrix is 

*L„(0 - X«.(0 + g / SVin,, S.) . (Al) 

where the plane wave is defined as 

X«,5.(r) = e^'^'-ns., (A2) 

further 

e{K) = K^/Mn , (A3) 
and K (M„) is the relative nn momentum (neutron mass). 

We expand the scattered wave as 

^^5.(^1 = 4vr E {fJ'S'i^'\j'\kJSu)^is',uisir)us'.', (A4) 

IJu'l'S' 

and the T-matrix as 

{k', S'u'\T\ii, Su) = J2 {^'^ l'S'u'\J^\k, ISu) T{!s',is{i^', k) ■ (A5) 
ii'j 

Then introducing the definitions 

^^hA^) = ^n,«j/i('^,0, (A6) 
TnK',K) ^ T/,. («:',«:), (A7) 



we obtain the equation, 

^zV(r) = / LuiM + I f[^"'''^"'"l, Tl,{.\.)\ . (A8) 



e[K) + if] — e[K' 
Let us write the function ^'^;j(r) as 

*?,aM = »"^^^. (A9) 
Then from Eq. (lASp . one obtains the following equation for the function wf, i{K,,r) 

'''' - 5v,i3i{Kr)+ / , , TijAK'.K). (AlO) 



' J e{K) + irj — e{K' 
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It can be shown that the new functions, defined as 



(All) 



satisfy the equations, 



K''^dK'jii{K'r) 



e{K) +iri - e{ 



^E^/',^«'^)f^n,.} • (A12) 



The unitary matrix Uf^ ^ is defined as 



COS e"' — sin e'^ 
sin e"^ cos e"^ 



(A13) 



where e"^ is the mixing angle. The relation between the T- and K-matrices reads 



I" 



I" a 



Here 



1 



(AM) 



(A15) 



Using Eq. ( ]A14I) in Eq ( ]A12[) . one obtains for the functions wi>^ja{n,r) the new set of equa- 
tions in terms of the K-matrix only 



KT 



It holds that 



cos 5i{jv{Kr)Ui',,^ + 



/I, a 



(A16) 
(A17) 



The functions w// ^{K,r) were computed from Eqs. (]A16|) by using the program phases 95 
Resolving Eq. (lAlip for the functions wf, i{K,r), one obtains 



W'A'^^'^) = E e'^^^il^/^a('^'0• 



(A18) 



Analogously, for the non-coupled channels we have 



*+(/€, r) 



I e 



KT 



(A19) 



where the function wi{K,r) satisfies the equation, 

wi{K,r) ^ i . , s f dn' 



cos (5; ijiinr) + 



e[K) — e K 



-K^ji{Kr)Ki{K) I . 



(A20) 
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Appendix B: The form factors arising after the Fourier transformation of the weak 
MECs due to the regulator ([2721) 



Here we present the form factors arising due to the regulator (12.21) after the Fourier trans- 
formation of the weak MECs containing one boson propagator with the mass m^. 



where 



TTXb 



inXB, A) - I^ixB, A) + XB ISixB, A) 

2 



= ^I2{xB,A) 

TTXb 



W2B = -^ vv^B 



Xb 



B 5 



/r(xB,A) = 

I^{XB,A) 

inxB,A) 

I2{XB,A) 



t sin (teg) ^(_t/A)2") 

1 + ^2 

' dt cos{txB) e("*/^^'") , 



1 + ^2 

: / dt sm{txB) e(-*/'^)") , 
Jo 

Xb = mBT , A = A/tub ■ 
For n=l, Eqs. (A.24)-(A.28) of Section A. 2 of Ref. [li] are reproduced. 



(Bl) 

(B2) 
(B3) 
(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
(B9) 
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